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In this work we calculate neutron and proton energy gaps in neutron star matter, using the Bonn 
meson-exchange interactions and a model-space approach to the gap equation. This approach allows 
a consistent calculation of energy gaps and single particle energies with the model-space Brueckner- 
Hartree-Fock (MBHF) method, without double counting of two-particle correlations. Neutron en- 
ergy gaps are calculated at zero and finite temperature. Proton energy gaps are calculated at beta 
equilibrium, and it is shown that the inclusion of muons has a significant effect. The results are com- 
pared with those of other works, and the implications for neutron star physics are briefly discussed. 



Superfluidity in nuclear matter has over the last 30 years received a great deal of attention, partly due to the 
possible impact on neutron star physics [1]. The existence of superfluid neutrons in the inner crust, and superfluid 
neutrons together with superconducting protons in the quantum liquid region of neutron stars must be considered as 
well established. Among the physical consequences of this qualitative picture, the most interesting is probably the 
role played by hadron superfluidity in cooling processes and in explaining pulsar glitches [2]. However, to assess the 
quantitative impact of superfluidity, precise knowledge of the size of the pairing energy gap is needed. Obtaining 



accurate results for this quantity in neutron star matter has proved to be a difficult problem. The most 



overwhelming difficulty is probably the consistent inclusion of medium effects in the effective pairing interaction, but 
even at the two-body level the results obtained so far are disturbingly scattered. This work is an attempt to clarify 
the situation at the two-body level, following an approach to the BCS gap equation introduced by Anderson and 
Morel [3] and applied to pairing in nuclear matter by Baldo et al. [4]. This approach to the gap equation allows us 
to calculate single-particle energies in the model-space Brueckner-Hartree-Fock (MBHF) approach in the same 
model space in which the gap equation is solved. In our calculations, we employ the meson-exchange potential 
models of the Bonn group [5], where the only parameters entering are physically motivated ones, like coupling 
constants and vertex cutoffs. In applications to finite nuclei and nuclear matter [6,7], the Bonn potentials have given 

satisfactory results. 

The structure of this work is as follows: In section 2 we reformulate the gap equation in terms of two coupled 
equations. We show how this reformulation allows us to renormalize the bare nucleon-nucleon (NN) interaction 
without introducing double counting of two-particle correlations in the gap equation. A brief description of the 
MBHF method is given, before we discuss the numerical solution of the gap equation. In section 3 we apply the 
method to pure neutron matter at zero and finite temperature, using the Bonn potentials of ref . [5]. In section 4 we 
investigate neutron star matter at [3 equilibrium, taking into account the appearance of muons at higher densities. 
In sections 5 and 6 we compare our results to those of similar works, and discuss the influence of the choice of 
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I. INTRODUCTION 



theoretical framework on the results. 



II. THEORY 



When the attraction in the partial wave is dominant, BCS theory [8,9] predicts a transition to a 
superfluid(neutrons)/superconducting(protons) state if a non-zero solution of the gap equation 
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can be found, where is the so-called gap function. Here k, k' are single-particle momenta^, Vk^k' is the 
component of the nucleon-nucleon interaction in momentum space, which is independent of the orientation of k and 
k', and = \J (e^ — tUpY + '^^ ^^'^ quasiparticle energy, where and t^p are the single-particle energies at k 
and kp respectively. The so-called decoupling approximation is assumed, that is, we take the Fermi surface to be 
sharp even in the presence of pairing correlations. This was found in ref. [10] to be a good approximation. Due to 
the decoupling approximation, the determination of single-particle energies is left as a separate problem which we 

will discuss later in this paper. 

The solution of equation (1) is complicated by the repulsive core of the nucleon-nucleon (NN) interaction, and a 
method for overcoming this difficulty has to be devised. One particular scheme is to replace the uncorrelated 
plane-wave basis used in the BCS approach with a basis in which short-range correlations are included, the 
so-called Correlated Basis Functions (CBF) approach, see Krotcheck and Clark [11] for an overview. Other authors 
have replaced the pairing matrix elements 14^' with the corresponding G-matrix elements [12,13], or with a pairing 
interaction constructed from Skyrme forces [14]. Another alternative was proposed in the work of Baldo et al. [4], 
where it was shown that the gap equation could be conveniently reformulated so that its role in accounting for 
short-range correlations becomes more transparent and the numerical treatment easier. As their description was 

rather brief, we will re-derive their result here. 
We split the summation in the gap equation (1) into two parts by introducing an arbitrary cutoff kc in momentum 

space: 

A, = - J2 V,,,^- J2 V,,,^. (2) 

k'<kc k'ykc 

Equation (2) is now taken as basis for an iterative procedure. As a first approximation, A'-'^-', we neglect the second 

sum in (2) and take 



k ^kk 2^ 

k'<kc 



Inserting this in the second sum in (2) gives a new approximation to the energy gap, A'-"'^-' 



k'<kc k'ykc k"<kc 

k'<kc k'<kc k">kc 

k'<kc k">kc 



We proceed by inserting the expression for A'-"'^-' in the second sum in (2) to obtain a new approximation A'-^-'. As 
long as Afe in the first sum is taken to be the exact solution of the gap equation for k < kc , ve need only iterate in 
the second sum. Continued iteration of this procedure gives an exact expression for the energy gap: 

Afe = — ^ {Vkk' — E ^kk"-^ — Vk"k' 
k'<kc k">kc 

k">kc k"'>kc K K K 

By defining 

Vkk' = Vkk' — Vkk"7rr< — Vk"k' (3) 

Zrjk" 

k">kc 



E E ^'='='"97^ — Vk"'k"T~^ — Vk"k' 



k">kc k"'>kc 
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^Throughout this work we set fi 
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the gap equation becomes 
A. = -^V>..,^. (4) 



k'<kc 



Furthermore, we he 



Vkk' — Vkk' — ^kk"-^ — {Vk"k' — Vk"k"'^rp, — 

k">kc k"'>kc 



Vk- 



Vkk' — ^kk"-^ — Vk"k' (5) 



k">kc 



This shows that an equivalent formulation of the gap equation (1) is given by the coupled integral equations 

Vkk' = Vkk' — ^ Vkk"-^ — Vk"k' (6) 



zEk' 

k">kc 

Ak = - J2 Vkk'^Ak', (7) 

k'<kc ^' 

which have to be solved simultaneously. We stress that there are no new approximations involved in these equations, 
and that they therefore represent a totally equivalent formulation of the gap equation (1). The role of the effective 
interaction V is to sum up all virtual transitions to momenta above kc- The advantage of this procedure is twofold. 
First, through the renormalization in equation (6) we obtain a more well-behaved pairing interaction, because the 
repulsive core of the bare interaction is integrated out. Secondly, we can explicitly see that double counting of 
two-particle correlations is avoided. Excitations to intermediate states above kc are included in V , while excitations 
to states below kc are included in the gap equation (7), see fig.l. Double counting will arise if the bare interaction 
V is replaced with any sort of effective interaction, e.g. the Brueckner G-matrix, which contains excitations to 
states below kc- This corresponds to taking kc = kp in equation (6), while at the same time taking kc = oo in 
equation (7), which is clearly inconsistent, see fig. 2. Thus caution is needed if pairing matrix elements are obtained 
from some sort of medium-renormalized interaction. One should then first find out which two-body correlations are 
included in the effective interaction, and then reformulate the gap equation so that the two-body correlations 

included in the effective interaction are excluded from the gap equation. 
Our calculation of pairing gaps is a two-step process: First we solve self-consistently the Brueckner-Hartree-Fock 
(BHF) equations for the single-particle energies, using a G-matrix defined through the Bethe-Brueckner-Goldstone 

equation 

G=V+V — %-G. (8) 

Here V is the NN interaction, Q is the Pauli operator which prevents scattering into intermediate states prohibited 
by the Pauli principle, Hq is the unperturbed Hamiltonian acting on the intermediate states and uj is the so-called 
starting energy, the unperturbed energy of the interacting nucleons. Methods to solve this equation are reviewed in 
ref. [6]. The single-particle energies for states ki (i encompasses all relevant quantum numbers like momentum, 
isospin projection, angular momentum, spin etc.) in nuclear matter are assumed to have the simple quadratic form 

^k. = ^+S., (9) 
2m* 

where m* is the effective mass. The terms m* and S, the latter being an effective single-particle potential related to 
the G-matrix, are obtained through the self-consistent BHF-procedure. The so-called model-space BHF-method is 
used [6,16], with a cutoff fcjvf > in momentum space. In this approach the single-particle spectrum is defined by 

ffe. = 77^ + «i (10) 
2m 

where m is the nucleon mass, and the single-particle potential Ui is given by 
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Ui = ^ {kikh\G{ui = eu, + ek^)\kikh) as , ki < kM 

kh<kF 

= 0, ki > kM, (11) 

where the subscript AS denotes antisymmetrized matrix elements. This prescription reduces the discontinuity in the 
single-particle spectrum as compared with the standard BHF choice k^j = kp. The self-consistency scheme consists 
in choosing adequate initial values for the effective mass and S. The G-matrix obtained is then used to calculate the 
single-particle potential Ui, from which we obtain new values for m* and S. This procedure continues until these 
parameters vary little. The nucleon-nucleon potentials are defined by the meson-exchange potential models of the 
Bonn group, in table A. 2 of ref . [5], versions A, B and C. By taking kc = in equations (6) and (7) we can 
calculate single-particle energies and solve the gap equation in the same model space. Although BCS theory within 
the decoupling approximation does not dictate the choice of single-particle energies, this model-space approach is a 
consistent and very satisfying method where all quantities of interest are calculated in a systematic manner within a 
model space containing the physically relevant degrees of freedom. We should point out that in the energy gap 
calculations we have used the single-particle spectrum defined by eqs. (10) and (11), and not the effective mass 

approximation (9). 

Having obtained single-particle energies, we proceed to solve eqs. (6) and (7). We choose kc = ^m, and solve the 
gap equation in the model-space in which the single-particle energies have been calculated. The solution of 
equations (6) and (7) can be simplified by noting that for k y kp, the dominating contribution to the quasiparticle 
energy Ek comes from the term (e^ — ^kp Y ■ This is clearly seen from fig. 3, where the quasiparticle energies are 
plotted for kp = 0.80 fm~^. Because is of the order 1 MeV, while Ek is of the order 100 MeV, we can safely 
neglect for k > k^j > kp. Thus equation (6) is decoupled from equation (7), and we can solve the linear equation 
for V by matrix inversion before proceeding to solve the gap equation. The gap equation in the continuum limit 

becomes 

Afe = — / rffcWfefe,-^. (12) 

Having obtained the effective interaction V , we solve equation (12) by iteration. As shown in fig. 4 the integrand in 
the gap equation has a relatively sharp peak at the Fermi momentum, so in the numerical treatment we have used 
50 Gauss-Legendre integration points in each of the intervals [0, kp + rj] fm~^ and [kp + rj, k^j] fm~^ in order to get 
numerically stable results, where rj = 0.1 fm~^ was found to be an adequate choice. Fig. 4 shows that it is necessary 

to use a rather large cutoff in the momentum space integration in the equation for V . 
To summarize: First we calculate renormalized single-particle energies using the MBHF method within a model 
space which size is determined by the parameter k^j- The next step is to evaluate the effective interaction V, which 
sums up pairing excitations to states outside the model space. Finally, the gap equation is solved with this effective 
interaction as pairing interaction within the model space. This method is both numerically efficient and stable, and 
also allows us to construct an effective pairing interaction without double counting of two-particle correlations. 

III. PURE NEUTRON MATTER 

As a first application we consider pure neutron matter at the relatively low densities found in the inner crusts of 
neutron stars. For simplicity we neglect the presence of finite nuclei in the crust, a presence which in refs. [18,19] 
was found to reduce considerably the neutron energy gaps compared with the pure neutron matter case. The 
calculation of energy gaps for realistic crust models is difficult, as one has to treat consistently the finite nuclei and 
the surrounding neutron gas. The main problem is to calculate pairing gaps in finite nuclei, and the validity of the 
local density approximations used so far has been questioned. In this section we therefore choose to focus on the 

more transparent case of pure neutron matter. 

A. Choice of model— space 

An important question is the choice of the parameter k^j which determines the size of the model space. Previous 
investigations for the MBHF-method [20] have shown that the results are stable in a region around k^j = 3.0 fm~^. 
Baldo et al. [4] found their results to be stable with respect to variations of kM, however, they did not give results 
for the ^M^dependence of the energy gap A^? = Ak=kF- In Hg- 5 we have plotted the energy gap at the Fermi 
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surface, Ap, calculated with the Bonn A potential, as a function of fcjvf for various choices of kp. The figure shows 
that the energy gap is virtually independent of fcjvf for in the interval 1.0 fm~^ < fcjvf < 4.0 fm~^. As fcjvf kp 
we get trouble due to our approximation AkykM — causing the energy denominator in equation (6) to almost 
vanish at kp. However, this problem is avoided in a fully self-consistent approach to equations (6) and (7), where 
Ai? gives a nonzero contribution to the energy denominator &i kp. For fcjvf 4.0 fm~^ large contributions from the 
repulsive core of V are included directly in the gap equation, and ultimately as fcjvf oo, the effective interaction 
becomes equal to the bare interaction. We notice that the numerical instabilities grow with increasing density. This 
is due to the increasing importance of the repulsive core. The results in this paper have been obtained with 
kM = 3.0 fm~^. Figures 6 and 7 show the ^M^dependence of the energy gap for three selected values of kp for Bonn 
B and Bonn C. In these calculations we included values of kM close to kp. It is seen that as k^j approaches kp the 
value of the energy gap decreases rapidly, due to the discontinuity in the single-particle spectrum. We see that the 

choice kjM = 3.0 fm~^ is right in the center of the stable region. 
The choice of upper limit for the numerical integration in equation (6) has also been shown to be important, and 
should according to ref. [17] be chosen greater than 50-60 fm~^ for the Paris potential, due to the complicated 
off-diagonal behaviour of this potential in momentum space. Accordingly, we have made sure that all results in this 

paper are stable with respect to the choice of upper limit for the integration in eq. (6). 



B. Energy gaps in pure neutron matter 



We have solved the gap equation in pure neutron matter with the potential models A, B and C of table A. 2 of ref. 
[5]. The results in fig. 8 show that the gap at the Fermi surface, Ap = A^^, is more or less the same for the three 
potential models. There is a slight difference in the maximum value of the energy gap, which is 
Afe^=o.80 = 2.71 MeV for Bonn A, Afe^=o.8i = 2.76 MeV for Bonn B and Afe^=o.8i = 2.75 MeV for Bonn C. The 
energy gap disappears for kp = 1.43 fm~^ with all potential models. These results should be compared with those of 
Baldo et al. [4], where a separable form of the Paris interaction gave the results Akp=o^sr^ = 2.76 MeV, A^^^i 4 = 0, 
and similar results were found with the Argonne V14 interaction. Our results agree with these results, as was to be 
expected, since all realistic nucleon-nucleon interactions reproduce the scattering data. Moreover, in this 
channel, only central force components contribute to the NN potential, so that more complicated off-shell effects 
from tensor forces are non-present. Thus at the two-body level all realistic NN-interactions should give similar 
results for the energy gap in pure neutron matter. The differences, apart from different choices for the single 
particle energies, may be traced to the off-shell behaviour of the interactions, which has been shown in ref. [4] to 
infiuence the value of the energy gap at kp. The advantage of the approach used here over that in [4] lies in the 
choice of MBHF instead of BHF single-particle energies, by which we avoid a discontinuity in the single-particle 

spectrum in the space where the gap equation (7) is solved. 
We are also interested in the gap function A^. To obtain the gap function for values of k above k^j we need to solve 
equations (6) and (7) self-consistently. This will also provide us with a check of the accuracy of the approximation 
Afe = 0, > kjM- The self-consistent solution of equations (6) and (7) is obtained by using the solution for 
Afe = 0, > kjM to calculate an initial approximation for A^^^kM from the gap equation (7). With this initial 
approximation we calculate a new effective interaction from equation (6) and then obtain a new approximation to 
Afe from equation (7). The procedure is continued until self-consistency is achieved. The results for the energy gap 
Ap were found to change insignificantly from the results with A^^^kM = 0, as seen from table 1. In figure 9 the gap 
function is shown for kp = 0.81 fm~^, calculated with the Bonn A potential. The gap function is seen to vary 
rapidly around kp and is negative over a wide range of momenta. However, contrary to the results with the Paris 
interaction in ref. [4], we find that the gap function decays rather quickly and shows little oscillation at large k. The 
difference is most likely due to the different off-diagonal momentum dependencies of the potentials. 

C. Correlated wave functions 

Another interesting quantity is the BCS pair wave function. For a given Fermi momentum kp, this function can be 

calculated from the gap function A^ using the expression 

/■OO A 

*^^^(r)=/ dkk''jo(kr)^ " (13) 
^ ' Jo ' ^ ' ^{ck - ck.r + Al' ^ ' 

which is the Fourier transform of the expectation value (c|a|j.^aljj.j^|c), with |c) denoting the usual BCS ground state, 
and ajj.^ creating a (normal) particle with momentum k and spin projection rris = +5 [4] The normalization 
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constant has been omitted. Fig. 10 shows this wave function together with the Brueckner two-particle wave 
function in the channel for two particles with opposite momenta at the Fermi surface, given by 

^io{r,kF) = joikpr) (14) 

^ Jkp ui - Ho(q,{)) 

where ui = -;^kp + 2S is the starting energy, Q{q, 0) is the Pauli operator, K is the center of mass momentum, and 
Go is the component of the G-matrix. Eq. (14) is a special case of the more general expression found in ref. [6]. 
We see that the behaviour of the two functions is very similar for small r, and this supports the view that the gap 
equation plays a part in accounting not only for the long range pairing correlations, but also for the short range 
two-particle correlations. A similar result was found by Baldo et al. [17] for pairing in the channel in symmetric 

nuclear matter. 



D. Finite temperature 



The gap equation at finite temperature T is 



A.(T) = -- r"rffcV2y,,,^^tanh f-f^) (15) 

where ks is Boltzmann's constant. To solve equation (15) we use the same approach as for T = 0. We make the 
simplification of ignoring the temperature dependence in V, because in the temperature range of interest 
ksT 1 MeV, while for k > kM , Ek is at least of order 100 MeV, and thus we can ignore the role of thermal 
excitations to states above kM- Furthermore, following common practice, we use "frozen" single-particle energies, 
that is, we calculate them within the MBHF approach at T = 0. According to refs. [21,22], this should be a 
reasonable approximation. We therefore calculate an effective interaction from equation (6) and use this to solve the 

gap equation (15) by iteration, just as for T = 0. 
In fig. 11 we show the temperature dependence of the energy gap at the Fermi surface, Ai?(T), calculated with the 
Bonn A potential. Also shown in fig. 11 are the critical temperatures estimated from the popular weak-coupling 

approximation (WCA) [23] 

fcsTc Pi 0.57Af(T = 0) (16) 

where Tc is the critical temperature. Our results indicate that this approximation is reasonably good, but that it 
systematically underestimates the critical temperature. The WCA seems to decrease in quality as the density 
increases, although not as drastically as the WCA for the energy gap and the condensation energy, which in ref. [4] 
were found to deviate considerably from the self-consistent results. That the WCA for the critical temperature 
turns out to be fairly accurate can be understood by considering its derivation, as found in e.g. ref. [23]. The main 
assumptions are that ksTc ^ £f and that it is sufficient to integrate over a small band of states near the Fermi 
surface. For neutron matter, the first assumption is reasonable, the last assumption is, however, dubious, as can be 

seen from fig. 4. However, Ap in eq. (16) is really the WCA for the energy gap, 

Af Pi 2e_Fexp ( 47r2-^^— ) (17) 

V Vkpkpkp J 

and if, as has been done here, it is replaced by the energy gap found by solving eq. (7), much of the error done in 

neglecting the integration over large momenta is compensated for. 

IV. PROTON PAIRING AT /3 EQUILIBRIUM 

We have also solved the gap equation for protons in neutron star matter at beta equilibrium, with and without 
muons. The proton fraction is determined by imposing the relevant equilibrium conditions on the processes 

e'+p^n + i^e (18) 
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e n +Vfj.+Ve (19) 
The conditions for [3 equilibrium require that 

IJ-n= IJ-P+ IJ-e (20) 

where jii is the chemical potential of particle species i, and that charge is conserved 

rip = Ue (21) 

where n,- is the particle number density for particle species i. If muons are present, the condition for charge 

conservation becomes 

rip = rig + (22) 
and chemical equilibrium in the process (19) gives the additional constraint 
IJ-e = IJ-iJ. (23) 

Throughout we have assumed that neutrinos escape freely from the neutron star. The proton and neutron chemical 
potentials are determined from the energy per baryon, calculated self-consistently in the MBHF approach. The 
electron chemical potential, and thereby the muon chemical potential, is then simply given by //g = A«n — The 

Fermi momentum of lepton type I = e, jj, is found from 

7 2 2 2 

Kp, = l^i - 

and we get the particle density using ri; = k'pj^tir'^. The proton fraction is then determined by the charge neutrality 
condition (22). Having calculated proton fractions and neutron and proton single-particle energies, we can solve the 
gap-equation for protons in the same manner as we did for neutrons in the previous section. In constructing the 
effective pairing interaction V for protons we neglect neutron-proton interactions, since the Fermi surfaces of 
neutrons and protons are widely separated. All results in this section have been obtained using the Bonn A potential. 
In fig. 12 the proton fractions obtained with and without muons are shown. The proton fraction is given as np/ns, 
where ns = n„ + np is the total baryonic number density. The muons enter the equation of state at a baryon 
density ns ^ 0.15 fm~^, and as a result the proton fraction increases. We also show the results of MBHF 
calculations of the proton effective mass in fig. 13. The faster decrease in the proton effective mass when muons are 
included is due to the larger proton density obtained with muons. Figure 14 shows the proton energy gap as 
function of the total baryonic density and in fig. 15 as a function of the proton Fermi momentum. In fig. 15 we also 
show the neutron energy gap calculated at /3 equilibrium. We see from these results that the inclusion of muons 
leads to a proton energy gap that decreases faster, but has a larger maximum value compared to the case with no 
muons. The faster decline of the energy gap is mainly due to the lower proton effective mass obtained with muons 
included, as shown in figure 13. The slightly larger maximum value of the energy gap is probably due to the increase 
in the proton density when muons are included, i.e. a proton "sees" more protons with which it can pair, and at the 
density where the proton gap reaches its maximum, this effect is greater than the effect of the increased density on 
the proton effective mass. The neutron energy gap is only slightly modified, due to the very low proton fraction at 
the baryon densities where neutron pairing is possible. We should mention that the presence of muons has little 
infiuence on the energy per particle and the resulting maximum masses for neutron stars, as discussed in [24]. 
To the best of our knowledge, there exists no previous calculations of the proton energy gap with muons included. 
For the results without muons, the most natural work to compare our results with is that of Baldo et al. [25]. There 
the Argonne V14 interaction [26] was used, and proton effective masses calculated in the BHF approximation. They 
found a maximum proton energy gap of approximately 0.8 MeV, and that \p[k^p) = at kp 0.9 fm~^. The 
differences between these values and those found in this work are probably due to the smaller proton effective mass 
ratios and smaller proton fractions obtained in ref. [25]. Amundsen and 0stgaard [10] obtained a maximum energy 
gap of 0.3 MeV using proton effective masses taken from Takatsuka [27], where a maximum proton gap of 0.5 MeV 
for kp = 0.5 fm~^ and a closing density of kp = 0.8 fm~^ was found. These low values are probably due to proton 
effective mass ratios in the range 0.5-0.7 being used in these two works. Chen et al. [15] obtained a maximum value 
of about 1.1 MeV within a variational (CBF) approach. However, as pointed out in that work, polarization effects 
will probably have a profound infiuence on the magnitude of the proton energy gap. Recent studies within the 
polarization-potential approach [28] indicate that medium effects can reduce the size of the energy gap by as much 

as a factor of three. 
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V. DISCUSSION 

The energy gaps calculated in this work are lower than those obtained with variational [10,15] and G-matrix [12] 
approaches. However, the results are not directly comparable, not only due to different choices of NN potentials, 
single-particle spectra, different calculations of proton fractions etc., but most importantly because they are 
obtained within different theoretical frameworks. As far as approaches using the Brueckner G-matrix is concerned, 
it is hard to see how these procedures can be justified in a strict manner. In perturbation theory, the G-matrix 
emerges as a natural reformulation of the original divergent perturbation series, but no similar reformulation of the 
pairing problem has been attempted within standard BCS theory. The CBF method used in refs. [11,15] employs 
basis states designed to handle short-range correlations, and at the pure variational level, the energy gaps obtained 
in this framework are larger than the bare interaction BCS results by some 0.5 MeV. It does not seem to be clear 
that double counting of two-particle correlations is avoided within the pure variational CBF approach, that is, that 
there is a clear division between the role of the pair correlation function and the correlation factors introduced in the 
CBF method. If medium polarization effects are omitted. Green's functions methods strongly indicate [29] that the 
pairing matrix elements of the bare nucleon-nucleon interaction should be used in the gap equation. The 
reformulation of the gap equation through the coupled equations (6) and (7) together with the wave functions in fig. 
10 provide evidence for the role played by the gap equation in taking care of the short-range correlations. The 
advantage of the model-space approach used in this work is that we can explicitly show that double counting is 
avoided. Furthermore, the gap equation is solved in a space where we use renormalized single-particle energies; free 
particle energies are only used in the construction of the effective interaction V . The renormalized interaction V 
makes the numerical treatment of the gap equation considerably easier, because the highly nonlinear part of the 
equation is separated from the almost linear part. We stress the need for extending the momentum space integration 
in the gap equation to large values of k. Thus the weak-coupling approximation used in ref. [13] is probably not 
appropriate, as also remarked in refs. [4,15], and tends to underestimate the energy gap. However, as they use a 
G-matrix or equivalent effective interaction in the construction of their quasiparticle interaction, this may 
compensate for the error caused by the weak-coupling approximation. The status of the results found in ref. [13] 

seems therefore to be uncertain. 

The inclusion of muons has a small, but significant effect on the proton energy gap. However, the proton energy 
gaps must be considered tentative, as we have neglected polarization effects. Furthermore, relativistic effects become 
important at the densities where proton pairing is possible. We hope to include both these effects in a future 
investigation. The proton fraction is of significant importance, and this again depends sensitively on the symmetry 
energy, which is a poorly known quantity. Albeit all these reservations, all calculations, including the one in this 
work, have so far been consistent with the existence of superconducting protons in neutron stars. 
Since all calculations, including the ones in this work, give neutron energy gaps of order 1 MeV, one can conclude 
that neutrons in the inner crusts of neutron stars are likely to form a superfiuid. Corrections to the simple and naive 
model of a uniform neutron gas interacting via two-body forces are probably not large enough to destroy neutron 
superfiuidity. However, as stated in the beginning of this paper, more realistic calculations are needed in order to 
study dynamical and thermodynamical effects in the inner crusts. The recent finding that the lattice nuclei in this 
region may take on non-spherical shapes [30,31] further complicates the picture. The impact of this fact on the 
specific heat of crust matter and the vortex pinning energy has been found to be significant [32]. The last quantity 

is of interest in models for pulsar glitches. 
The proton energy gaps obtained in this work have recently been used in a calculation of rates for the modified 
URCA processes [33]. Not surprisingly, the rates were found to be significantly suppressed by proton 
superconductivity. The analysis of Page [34] indicates that it is hard to distinguish fast from slow cooling scenarios; 
in both cases agreement with the observed temperature of e.g. Geminga is obtained if baryon pairing is present in 

most, if not all of the core of the star. 

VI. CONCLUSION 

The reformulation of the BCS gap equation advocated by Baldo et al. [4,17] allows for a consistent treatment of 
pairing correlations without double counting. One also gets the advantage of using a renormalized interaction, so 
that the main effects of the repulsive core of the NN interaction can be integrated out before proceeding to solve the 
gap equation. The additional consistency in solving the gap equation in the same model space in which the 
single-particle energies are calculated is also an advantage of this method. However, medium polarization effects on 
the effective pairing interaction must be taken into account, and work in this direction is in progress. From our 
point of view, the model-space approach to pairing correlations seems to offer possibilities also for a consistent 
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treatment of medium polarization through the inclusion of particle-hole ring diagrams. 
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One of us, (MHJ), thanks the NFR and the Istituto Trentino di Cultura, Italy, for financial support. 



[1] S. L. Shapiro and S. A. Teukolsky, Black Holes, White Dwarfs, and Neutron Stars, The Physics of Compact Objects (Wiley, 

New York, 1983) p. 290 
[2] G. Baym, C. J. Pethick and D. Pines, Nature 224 (1969) 673 
[3] P. W. Anderson and P. Morel, Phys. Rev. 123(1961) 1911 

[4] M. Baldo, J. Cugnon, A. Lejeune and U. Lombardo, Nucl. Phys. A515 (1990) 407 
[5] R. Machleidt, Adv. Nucl. Phys. 19(1989) 189 

[6] M. Hjorth-Jensen, T. T. S. Kuo, and E. Osnes, Phys. Rep. 261 (1995) 1 

[7] R. Brockmann and R. Machleidt, Phys. Rev. C 42 (1990) 1965 

[8] J. Bardeen, L. N. Cooper and J. R. Schrieffer, Phys. Rev. 108 (1957) 1175 

[9] L. N. Cooper, R. L. Mills and A. M. Sessler, Phys. Rev. 114 (1959) 1377 
[10] L. Amundsen and E. 0stgaard, Nucl. Phys. A437(1985) 487 
[11] E. Krotscheck and J. W. Clark, Nucl. Phys. A333 (1980) 77 
[12] E. 0stgaard, Z. Physik 243(1971) 79 

[13] T. L. Ainsworth, J. Wambach and D. Pines, Phys. Lett. B222 (1989) 173 
[14] M. E. Siang and T. T. S. Kuo, Nucl. Phys. A481 (1988) 294 

[15] J. M. C. Chen, J. W. Clark, R. D. Dave and V. V. Khodel, Nucl. Phys. A555(1993) 59 
[16] T. T. S. Kuo and Z. Y. Ma, Phys. Lett. B127 (1983) 123 

[17] M. Baldo, U. Lombardo, I. Bombaci and P. Schuck, Phys. Rep. 242 (1994) 159 

[18] R. A. Broglia, E. V. De Blasio, G. Lazzari, M. C. Lazzari and P. M. Pizzochero, Phys. Rev. D 50 (1994) 4781 

[19] 0. Elgar0y, L. Engvik, E. Osnes, E. V. De Blasio, M. Hjorth-Jensen and G. Lazzari, LANL preprint xxx.nucl-th/9509035 

[20] H. Q. Song and T. T. S. Kuo, Phys. Rev. C 43 (1991) 2883 

[21] A. Lejeune, P. Grange, M. Martzolff and J. Cugnon, Nucl. Phys. A453 (1986) 189 

[22] M. Baldo, I. Bombaci, L. S. Eerreira, G. Giansiracusa and U. Lombardo, Phys. Lett. B215 (1988) 19 

[23] E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, Part 2 (Pergamon, Oxford, 1980) 

[24] L. Engvik, G. Bao, M. Hjorth-Jensen, E. Osnes, and E. Ostgaard, submitted to Ap. J. 

[25] M. Baldo, J. Cugnon, A. Leujeune and U. Lombardo, Nucl. Phys. A536 (1992) 349 

[26] I. E. Lagaris and V. R. Pandharipande, Nucl. Phys. A359 (1981) 331 

[27] T. Takatsuka, Prog. Theor. Phys. 50 (1973) 1754 

[28] J. Wambach, T. L. Ainsworth and D. Pines, in Neutron stars; Theory and observation, eds. J. Ventura and D. Pines 
(Kluwer, Dordrecht, 1991) 

[29] A. B. Migdal, Theory of finite Eermi systems and applications to atomic nuclei (Wiley, New York, 1967) 

[30] C. P. Lorenz, D. G. Ravenhall and C. J. Pethick, Phys. Rev. Lett. 70 (1993) 379 

[31] K. Oyamatsu, Nucl. Phys. A561 (1993) 431 

[32] G. Lazzari, E. V. De Blasio, Z. Phys. A353 (1995) 13; 235 

[33] 0. Elgar0y, L. Engvik, E. Osnes, E. V. De Blasio, M. Hjorth-Jensen and G. Lazzari, submitted to Phys. Rev. Lett. 
[34] D. Page, Ap. J. 428 (1994) 250 



10 



a) 



b) 





\ 1 















































— 

/ 




\ 1 


. — 

\ 



























































FIG. 1. Two-particle excitations involved in pairing. Those in a) are included in the gap equation (6), and those in b) are 
included in the effective interaction (5). 
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FIG. 2. If we take fc^ = A;_p in the effective interaction, while fee = oo in the gap equation, we get the situation shown above, 
where the two-particle excitation in a) and b) is included in both equation (5) and in equation (6). 
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FIG. 3. Quasiparticle energies in pure neutron matter at kp = 0.80 fm~^ calculated with the Bonn A potential and 
km = 3.0 fm~^ . 
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FIG. 4. The integrand in the gap equation for kp = 0.80 fm and k = 0.90 fm , calculated with the Bonn A potential. 
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FIG. 5. Energy gap A_p with the Bonn A potential plotted for various values of /cm- The numbers next to the lines are the 
corresponding values of kp. 
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FIG. 6. Dependence upon /cm of the energy gap calculated with Bonn B for kp = 0.50, 0.80 and 1.10 fm 
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FIG. 7. Dependence upon /cm of the energy gap calculated with Bonn C for kp 
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FIG. 8. The energy gap A_p calculated with the Bonn A, B and C potentials. 
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FIG. 9. Gap function A(k) = Ak calculated with the Bonn A potential for kp = 0.81 fm~^. 
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FIG. 10. The wave functions and *o for kp = 1.1 fm"^ 
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FIG. 11. Temperature dependence of A_p for the values of A;_p (in fm~^) indicated in the figure. The corresponding 
weak-coupling estimates for the critical temperatures are also indicated. 
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FIG. 14. Proton energy gap at /3 equilibrium, shown as a function of baryon density. 
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FIG. 15. Proton and neutron energy gaps as functions of the respective Fermi momenta. 
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0.50 
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1.882 

0.70 
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2.714 

2.711 

0.85 

2.693 

2.689 

1.00 

2.304 

2.301 

1.10 

1.800 

1.798 

1.30 

0.522 

0.522 

TABLE I. Energy gaps calculated with Bonn A. Here A^' has been obtained with the approximation Akyk^ = 0, while 

(2) 

Ap is the complete self-consistent solution. 
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